Universal thermal conductivity in the vortex state of cuprate superconductors 
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We formulate an effective low energy theory for the fermionic excitations in d-wave superconductors 
in the presence of periodic vortex lattices. These can be modeled by an effective free Dirac Hamil- 
tonian with renormalized velocities and possibly a small mass term. In the presence of random 
nonmagnetic impurities this will result in universal (i.e. field and disorder strength independent) 
thermal and spin conductivities with values different from those occurring in the Meissner state. 



At low energies physics of the cuprate superconduc- 
tors is dominated by fermionic excitations in the vicinity 
of the four nodes of the dx2^y2 superconducting order 
parameter. FormaUy, these low energy excitations can 
be described as four species of relativistic massless Dirac 
fermions. Perhaps the most spectacular manifestation of 
these Dirac-like excitations is the appearance of univer- 
sal conductivities at low temperatures. This phenomenon 
was predicted first by Lee |l| for electrical conductivity 
and later extended to spin and thermal conductivities 
P,B|. "Universal" in this context means that as a result 
of the linear dispersion of Dirac quasiparticles conductiv- 
ities become independent of the scattering rate below cer- 
tain temperature scale (which itself is non-universal and 
set by the scattering rate). Measurement of these uni- 
versal conductivities thus provides information about the 
intrinsic properties of the underlying clean system. Ther- 
mal conductivity has a special significance since its uni- 
versal value is unaffected by vertex and Fermi liquid cor- 
rections 0] . Experimental measurements of thermal con- 
ductivity at sub-Kelvin temperatures confirmed the exis- 
tence of the universal regime in samples of YBa2Cu306.9 
(YBCO) g and BiaSraCaCuaOs || and provided values 
of the Dirac anisotropy ratio aiy — vf/va in agreement 
with other probes, such as the angle-resolved photoemis- 
sion. To date, no observation of universal conductivity 
in charge or spin channel has been reported. 

In the present paper we argue that, under certain con- 
ditions, universal conductivities may also appear in the 
vortex state of a d-wave superconductor at low tempera- 
tures. In this case "universal" means independent of both 
the quasiparticle scattering rate and the applied mag- 
netic field. The physics of this result is based on the 
recent body of work [[7|-|l3[ that resulted in detailed un- 
derstanding of the quasiparticle dynamics in the presence 
of a vortex lattice. Within the simplest linearized Dirac 
model the quasiparticles are found to form Bloch bands 
in the periodic scalar and vector potentials created by 
the vortex lattice. One important conclusion is that the 
original Dirac node is preserved in the vortex state [QH] 
albeit with renormalized quasiparticle velocities vp and 
VA- Thus, at the lowest energies, one may describe the 
quasiparticle dynamics in the vortex state by an effec- 
tive free Dirac Hamiltonian with renormalized parame- 



ters. Since the T — > behavior is entirely determined by 
this low energy Hamiltonian, computation of conductiv- 
ities follows essentially the same path as in the Meissner 
state Kn and yields the aforementioned universal be- 
havior. 

The above sketch captures the essential physics of our 
result and would be exact within the linearized model. 
The reality is somewhat more complicated. Calculations 
beyond the linearized model |jiq,na,n3[ indicate that in- 
ternodal scattering and nonlinear terms in the full Hamil- 
tonian can modify the band structures discussed above 
by displacing or creating additional Dirac nodes or pro- 
ducing small gaps at the Dirac point. Such "massive" 
Dirac fermions have been argued to give rise to quan- 
tized thermal Hall conductivity n^y |12,M. We demon- 
strate below that, surprisingly, even such massive Dirac 
fermions give rise to universal longitudinal conductiYities, 
provided that the scattering is in the unitary limit. 

Thermal conductivity in the vortex state of a d-wave 
superconductor has been addressed previously by number 
of authors ||l4|-[l7| . In these works the effect of vortices 
has been treated within semiclassical "Volovik" approxi- 
mation [Q which tends to capture the essential physics 
in many situations |l9-pl|] but whose domain of valid- 
ity remains under debate Q. Here, for the first time, 
we present a fully quantum treatment of the longitudi- 
nal thermal and spin conductivities in the vortex lattice. 
Such quantum treatment is essential when addressing the 
behavior of very clean systems in the T ^ limit. 

Our starting point is the BogoHubov-deGennes (BdG) 
Hamiltonian linearized near a single node of the d-wave 
order parameter, 

H == vf{Px + ax)cr3. + VAiPy + ay)ai 

+ VpVsx + U{v)(T3. (1) 

Here Ci are the Pauli matrices, p = — ^V, and we have 
already performed the singular gauge transformation M] 
to "unwind" the phase of the superconducting order pa- 
rameter A = Aqc"*'''. Dirac fermions are coupled to 
the "Berry" gauge field a = i(V(/?A — V(/3_b), and to 
the "Doppler" gauge field Vg = \{'^'fA + V(p_B - 2eA/c) 
which formally enter Hamiltonian (Q) as vector and scalar 
potentials respectively. In a static vortex lattice these 
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FIG. 1. Low energy band structure for Dirac quasiparticles 
with an = 4 in the presence of square vortex lattice (solid 
lines) and in zero field (dashed lines) . Lower panel: associated 
low energy density of states, luh ~ 'k^Jvfv/\{H/2^o). 



two gauge fields are determined through gradients of the 
phase fields LpA and lpb- These represent the singular 
parts of the phase coming from two sublattices A and B 
of the vortex lattice, and satisfy ip = ipA + <^b- The last 
term in Eq. (|l|) represents the effects of nonmagnetic im- 
purities. Since disorder couples to the electron density, 
it is invariant under the singular gauge transformation. 
In the absence of disorder and for periodic vortex lat- 
tices Hamiltonian (|l]) can be diagonalized using stan- 
dard band structure techniques [0J|,|ll|. As illustrated 
in Fig. |l|, for inversion symmetric vortex lattice, spec- 
trum remains Dirac-like at the lowest energies; the main 
effect of the vortex lattice is to renormalize the Dirac 
velocities vp and va ■ Once the presence of the nodal 
points is established H,@,n2|, it is straightforward to 
understand the linear dispersion. The wavefunctions 
can be written in the Bloch form, 'l'k(r) ~ e''^'"xk(r), 
where Xk(r) is a two component spinor periodic on the 
unit cell. Xk(r) is an eigenstate of the Bloch Hamilto- 
nian 7i(k) = e^^^'^Hoe^^'^ , where Ho is Hamiltonian 
(|l|) with the disorder term set to zero. At zero energy 
there is a doublet of degenerate states, |xo )' satisfying 
H(0)|xo ) = 0. These zero energy states provide a con- 
venient basis for a perturbative expansion in the vicinity 
of the node: we may write 



7i;(k) == 7i:(0) + {vFk^<J3 + VAkydi), 



(2) 



and treat the second term as a small perturbation to 7i(0) 
|l2| . We find that, near the nodal point, a simple first or- 
der degenerate perturbation theory yields very accurate 
results. For the energy we obtain 



E, 



(1) 



± 



\j{vFkx)^ + {vi\kyY 



(3) 



where vp and wa are renormalized Dirac velocities to 
be specified shortly. Higher order terms produce O(k^) 
corrections to the energy and are unimportant at small 
k, in agreement with Fig. |l|. 

The above considerations imply that very accurate de- 
scription of the low energy excitations of the system in 
the absence of disorder can be achieved by projecting its 
Hamiltonian (g) onto the subspace spanned by its two 
zero energy eigenstates. We now include the effects of 
disorder by projecting the full Hamiltonian (P onto this 
subspace. Formally this is accomplished by introducing a 
projector V = Ek,.=± IM(kj^l with |k,±) = e^'^-'lxo )■ 
The effective low energy Hamiltonian 7i°* = VHV reads 

^kk' = f^kk' [vFkx{ni ■ t) + VAky{n2 ■ t)] + Ukk' • r, 

where Ti are Pauli matrices acting in the (/x, v) space 
and ni.2 are constant vectors defined by [rii • t]^j, = 
(k/i|(73|ki/) and [n2 ■ t]^^ — (k/i|(Ti|kj/). Similarly Ukk' 
is defined by [Ukk' ■''']tJ.i' — (k/i|cr3t/(r)|kV). Making the 
usual assumption C/(r) = uq J2a <^(i"~i"q) where {Va} are 
random impurity positions, we obtain 
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U=uoJ2^'"'"'^'''''''''KAG), 



(4) 



(x'^ 



oPse" 



ir G 



Xo) 



Summation over the 



withZ^^^(G) 

reciprocal lattice vectors G arises because of the peri- 
odicity of Xo(r) and will in general complicate disorder 
averaging. However, numerical evaluation ofU^^{G) in- 
dicates that the sum in Eq. (0) is dominated by the uni- 
form (G = 0) term. In the following we therefore drop 
all but this imiform term, i.e. U^^{G) -^ (5G=o[ni • t]^^. 
This approximation allows for simple disorder averaging. 
For vortex lattices with inversion symmetry it is easy to 
show |12] that ni • n2 = 0. In such a case, in the absence 
of disorder, 7i°^ has spectrum (|^) with vp —VF\ni\ and 
VA — VA\n.2\- It is then also possible to choose as a 
basis such linear combination of the degenerate states 
IXo ) that vpi'n.i ■ r) = vpT^ and va{^2 ■ t) = vati and 
write the effective low energy Hamiltonian as 



'Ikk' 



5kk' [vFkxT3 + VAkyTi + mDT2] + C/kk'T3. (5) 



Here t/kk' = Mo|ni| X]a ^"^"^"^ "^ ^ ^^^ we added by 
hand the mo "mass" term to model the small gap which 
according to Refs. |10||l^ can open up at the Dirac point 
as a result of internodal scattering or nonlinear terms ne- 
glected in iM). In the absence of disorder the spectrum of 
Hamiltonian (||) is E^ = zty^iypk^)^ + {vAkyY + mjy. 

The effective Hamiltonian (|^) is valid at energies E <C 
u!h = IT ^J V pv a{E[ / 'i.^o) where $o — hc/2e is the flux 
quantum. Taking the YBCO values (6| vf = 2.5x lO^m/s 
and ao = 14 we find u>h = 25K-\/i//lT: Hamiltonian 
(o) should be valid at sub-Kelvin temperatures relevant 
to the low-T heat conduction experiments P,|q,pi|| . 

Around the single Dirac node the bare Matsubara 
Green's function can be written as 



^o(k, iw„ 



veKts + VAkyTi + moT2 



2 ^ 

D) 



(6) 



Ns 



27r 



duj K [uj) (11) 



wherein — (2n+l)7rT and e^ — {vphxY + {v\kyY. The 
impurities alter the bare Green's function by introducing 
a Matsubara self-energy S(ia;„). In the spirit of Ref. Q] 
we assume that all but the scalar component of S(ja;„) 
can be neglected or absorbed into dispersion or pairing 
fE2i. Hence the dressed Green's function becomes 



where s = 1/2 is the coupling constant for spin current 
and 



K\lo) 



_fk_ 



({}|Tr[C7;4t(k,c.)r3Ct(k,c^)r3] 
-vlTT[g:UK^)rig:UKu;)n]), (12) 



g(\i,iLu„) = go{'k,iujn - S(iw„)). 



(7) 



Retarded Green's functions are obtained by analytically 
continuing 5i-et(k, w) = t/(k, ia;„ ^ lu + iS) and the im- 
purity scattering rate is defined as 7(0;) = — Im Ei-et(w). 
Within the self-consistent t-matrix approximation the 
self energy is given by p2|,p3[ 



E(iw„) = TgQ{iuJn)/[c'^ - gl{iujn)], 



(8) 



with r = ni/irpo, rii the impurity density, po the normal 
state DOS; c = cotJo with 5q the scattering phase shift, 
and go{iujn) = (27r/Oo)~^A /"y] v TrQ {]<., iun) , N the num- 
ber of Dirac nodes. Eq. (||) self-consistently determines 
the frequency dependent scattering rate 7(w)- At low 
temperatures we are interested in 70 = 7(0; —> 0). In the 
Born limit (c ^ 1) we find 



ll 



^2^^p2^-A.c-v^v^/MT 



(9) 



where A is the upper cutoff of the order of maximum 
superconducting gap. In the massless case this equation 
always has a real solution for 70 implying finite DOS as 
tj — > and universal conductivities, albeit below expo- 
nentially small temperatures 0]. When mjj > there 
is no real solution below the critical impurity concentra- 
tion n1 — 2T:'^c^VFVAPo/[Mhi{K/mD)\'- in the massive 
case weak disorder cannot fill in the gap and produce 
universal conductivities. 

In the unitary limit [c —f 0), we find 



7o 



^ ~ TT^VpV^r 



A/" In 



A2 



I'o 



'D 



(10) 



This equation has real solution for arbitrarily small impu- 
rity concentrations. In the presence of unitary scattering 
the system will exhibit finite DOS at lo = and, as we 
show below, universal conductivities even when mjj > 0. 
Physically this is a consequence of impurity bound states 
forming inside the gap p3[ . Overlap between these states 
leads to formation of impurity band capable of carrying 
the quasiparticle current. 

We proceed by computing the spin conductivity a^ 
which is simply related to the heat conductivity by 
Wiedemann-Franz law |,|l|l, k/T = {TT^k%/3s'^)a' . 
From now on we shall assume that 70 > 0, which ac- 
cording to the discussion above is guaranteed in the uni- 
tary scattering limit. Neglecting vertex corrections H] 
the static spin conductivity reads 



with g'^'^^{k,uj) = ImQj-ctiKuj). 
can take —dnF{Lo)/dLu — > S{uj). 



In the limit T 
Substituting 



^;ct(k,0) 



70 - ■mDiT2 



'D 



into Eq. (|l2|) and performing the traces we find 



K\0) 
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(e2-Hm|3+72) 
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we 



(13) 



(14) 



The remarkable feature of this result is that the Dirac 
mass and the scattering rate enter only in the combina- 
tion 77i|j -l-7o- This is somewhat counterintuitive since in 
the single particle spectrum the two tend to have oppo- 
site effects: mass depletes the low energy DOS while 70 
enhances it. 

By power counting the integral in Eq. m% is seen to be 
independent of m|, + 7q , implying universal conductivity 
just as in the massless case. We obtain 



Ms^ vl 



47r2 



VpVA 



and by Wiedemann-Franz law 
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(15) 



(16) 



= vp/vA S> 1 is the renormalized Dirac 
Direct computation of k yields identical re- 



where an 

anisotropy. 

suit. 

At fields Hci <C iif <C i?c2 the only scale in the problem 
is the intervortex separation Z, implying that to leading 
order an is field independent |24[| y et distinct from H = 
value. Thermal conductivity (|lq) is universal: it is in- 
dependent of disorder and magnetic field H. Numerical 
calculations within linearized model RSO] show that 
a_D > ao and thus we expect the quasiparticle contribu- 
tion to the thermal conductivity to be enhanced in the 
presence of vortices [p5| . 

Fig. y summarizes our result for the field dependence 
of the longitudinal thermal conductivity in the T —^ 
limit in the presence of periodic inversion symmetric 
vortex lattice. In the Meissner state {H < Hd) mag- 
netic field is excluded from the sample and thermal 
conductivity assumes its zero field universal value g] 
Ko/T = {kg/3){vF/vA)- Above Hd formation of the 
vortex lattice causes renormalization of the Dirac ve- 
locities and possible gapping or creation of additional 
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FIG. 2. Schematic plot of the universal quasiparticle ther- 
mal conductivity k in the vortex state at T = 0. The con- 
tribution from TV species of Dirac fermions is proportional 
to the renormalized Dirac anisotropy old > old S> 1. In the 
shaded region just above Hd the magnetic field is nonuniform 
(i 2> A_l) and our model does not apply. 



nodes. The associated universal conductivity is k/T = 
{N'k'^ /12){vf /va)- It is worth noting that this type of 
behavior has been recently observed in the ultra-pure sin- 
gle crystals of YBCO |6[. 

We conclude by discussing the range of validity of our 
considerations. The crucial assumption is that of per- 
fectly periodic Bravais vortex lattice with inversion sym- 
metry. Within the linearized model this guarantees the 
existence of Dirac quasiparticles at low energies [p|,p^. 
This assumption should be well satisfied in the ultra- 
pure single crystals of YBCO. Inclusion of internodal 
scattering and nonlinear effects |10| , p^ neglected in the 
linearized model, in general can produce small gaps at 
the Dirac points or introduce additional Dirac nodes in 
the spectrum [1^] . With the mass term present our effec- 
tive Hamiltonian (ra) is sufficiently general to include any 
such effects. A remarkable result is that, in the presence 
of unitary scatterers (or sufficiently high concentration of 
Born impurities) even massive Dirac fermions give rise to 
the universal conductivities. We note in passing that this 
implies universal conductivities for situations in which 
the Dirac nodes have been gapped by other mechanisms, 



such as in d^2_ 



idxy state. 



This conclusion holds 



within the self-consistent t-matrix approximation which 
corresponds to the saddle point analysis in the replicated 
field theory treatment of Rcf. |g] . Ultimately, at the low- 
est energies, fluctuations around this saddle point will 
drive the DOS to zero. It appears, however, that this 
regime has not yet been accessed experimentally, except 
perhaps in samples of underdoped YBCO [p7| . 

Linearized models indicate that for bare Dirac 
anisotropics a^ > 10 the band structure in the vortex 
lattice becomes essentially one dimensional |7|,|8|,|ll| , i.e. 
aD ^ cx). In such a case our effective Hamiltonian (0) 
would have vanishing domain of validity. Calculations us- 
ing the full BdG Hamiltonian ||lQ,tt3| however show that 
the renormalization of the anisotropy is much less severe. 

In less pure samples we expect disorder to destroy pe- 
riodic vortex lattice. It has been argued that universal 
field independent conductivity also emerges in the ran- 



dom vortex arrays at high fields H > Hc2{T/Tcf [|T|,||. 
This type of behavior was indeed observed at T > ?>K 
pO| . At sub-Kelvin temperatures, on the other hand, 
it was found in less pure samples that n/T ^ \/H [ pi[ , 
consistent with semiclassical treatments |lj,nq] which ne- 
glect quasiparticle scattering from vortices. From this 
perspective our result for clean samples sketched in Fig. 
represents a nontrivial universal limit with a simple pre- 
diction that is testable by experiment. Although the size 
of the jump at Hd depends on the details its presence is 
a robust qualitative prediction of the present theory. 

The authors are indebted to R. Hill, C. Lupien and L. 
Taillefer for discussing their data on ultra-pure YBCO 
crystals prior to publication, and are grateful to M.-R. 
Li, A. Melikyan and Z. Tesanovic for discussions. This 
work was supported in part by NSF grant DMR-9415549 
(O.V.) and by NSERC (M.F.) 
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